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Abstract 

We investigate the implication of the recent discovery of a Higgs-like particle in the first phase 
of the LHC Run 1 on the Inert Higgs Doublet Model (IHDM). The determination of the Higgs 
couplings to SM particles and its intrinsic properties will get improved during the new LHC Run 2 
starting this year. The new LHC Run 2 would also shade some light on the triple Higgs coupling. 
Such measurement is very important in order to establish the details of the electroweak symmetry 
breaking mechanism. Given the importance of the Higgs couplings both at the LHC and e'’"e“ 
Linear Collider machines, accurate theoretical predictions are required. We study the radiative 
corrections to the triple Higgs coupling hhh and to hZZ, hWW couplings in the context of the 
IHDM. By combining several theoretical and experimental constraints on parameter space, we 
show that extra particles might modify the triple Higgs coupling near threshold regions. Finally, 
we discuss the effect of these corrections on the double Higgs production signal at the e+e“ LC 
and show that they can be rather important. 
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I. INTRODUCTION 


The discovery of a new particle with a mass around 125-126 GeV in the search for the 
Standard Model (SM) Higgs boson [T] was announced simultaneously by the ATLAS and 
CMS collaborations in July 2012 [2H7]. Since then, more data has been taken and analyzed 
at the LHC. One of the primary goals of the Higgs groups at the LHC is now to study the 
properties of this new resonance and determine if it is indeed the state predicted by the SM. 

With this new discovery, a program of precision measurement involving the Higgs boson 
has just started and will get improved with the new run of LHC and future Linear 

Collider (LC). In fact, the 7 © 8 TeV data allow the measurement of the Higgs couplings to 
gauge bosons and r+r“ with about 20-30% of precision while the Higgs couplings to bb and 
ti still suffer large uncertainties of about 40-50%. All these measurements will be improved 
with the new run of LHC at 13-14 TeV and the future LC. 

In order to conhrm that the discovered Higgs-like particle is the SM Higgs responsible 
for the electroweak symmetry breaking, we need to know all its couplings to SM particles 
with accurate precision and also measure the trilinear and quartic self-couplings of the Higgs 
in order to be able to reconstruct the scalar potential. In this regards, the LHC with its 
high luminosity option have also the capability of measuring the SM triple Higgs couplings 
through one of the following channels gg ^ hh ^ bb'j'j,bbT~^T~,bbW^W^* [SI E]- The 
e'''e“ LC, which will provide some high precision measurement of the Higgs mass and its 
properties such as couplings to SM particles and quantum numbers, would also be able 
to perform SM triple Higgs coupling through Zhh (double Higgs-Strahlung) and 

e’''e“ —)■ Ve'i^ehh (WW fusion) with more than 700 GeV center of mass energy with better 
precision [10]. In the double Higgs-Strahlung process Zhh the Z boson will be reconstructed 
from l^l~ or qq pairs, while for the WW fusion process the two Higgs can be reconstructed 
from bhbh or bb and W~^W~. 

The discovery of this Higgs-like particle resonance opens a new era in elementary particle 
physics and leads to several theoretical and phenomenological studies on Higgs physics both 
in the SM and beyond. One of the very simplest extension of SM is the IHDM proposed, 
more than three decade ago, by Deshpande and Ma m for electroweak symmetry breaking 
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purpose. Recently, the IHDM model has been very attractive because it provides a dark 
matter candidate [HI US], generates tiny neutrino masses [Hj and also solve the naturalness 
problem [15]. Phenomenology of IHDM have been extensively studied during last decade 

[MiiiHiini. 

The aim of this paper is to study the effect of the one loop radiative corrections to the 
triple Higgs coupling hhh as well as hZZ coupling in the framework of the IHDM. We will 
also compute the well known SM radiative corrections to the triple Higgs coupling as check 
of our procedure. Once these effect are well studied, we then proceed to the evaluation of 
radiative corrections to the double Higgs Strahlung process e’''e“ —)■ Zhh. For this purpose, 
we will apply an on-shell renormalization scheme to evaluate these one-loop corrections. 
For the numerical evaluation, we will take into account all the theoretical and experimental 
constraints on the scalar sector of the Model. 

The paper is organized as follow: in the second section we introduce the IHDM model 
and describe the theoretical and experimental constraints that the model is subject to. 
In the third section we introduce the on-shell renormalization scheme for the triple Higgs 
coupling hhh and hZZ in the IHDM and present our numerical results in the fourth section. 
Numerical analysis of the double Higgs-strahlung is presented in the hfth section. Our 
conclusion is given in the last section. 

II. THE INERT HIGGS DOUBLET MODEL 

A. The Model 

The IHDM is one of the most simplest models for the scalar dark matter, a version of 
a two Higgs double model with an exact Z 2 symmetry. The SM scalar sector is extended 
by an inert scalar doublet H 2 which can provide a stable dark matter candidate. Under Z 2 
symmetry all the SM particles are even while H 2 is odd and it could mix with the SM-like 
Higgs doublet. We shall use the following parameterization of the two doublets : 

( 

Hi = \ , H2= \ 

\^^(v + h + i0>) ) + */!») 
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with and are the Nambu-Goldstone bosons absorbed by the longitudinal component 
of and respectively, v is the vacuum expectation value (VEV) of the SM Higgs 
Hi. Within the IHDM the scalar doublet H2 does not couple with the SM fermions and 
therefore the i^f 2 -fermions interaction are present only through mixing with Hi. The most 
general renormalizable, gauge invariant and CP invariant potential is given by : 

V = + /i2|hf2p + + A2|hf2|^ + A3|iPip|ip2p + A4|ipjip2p 

+ + h-c} (2) 

In the above potential there is no mixing terms like Hi 2 {HiH 2 + h.c) because of the unbroken 
Z2 symmetry. By hermicity of the potential, all Aj, f = 1, • ■ ■ , 4 parameters are real. The 
phase of A 5 can be absorbed by a suitable redehnition of the helds Hi and H2, therefore the 
scalar sector is CP conserving. After spontaneous symmetry breaking of SU{2)l 0 U{1)y 
down to U{1)q, the spectrum of this potential will have hve scalar particles: two CP even 
H^ and h which will be identihed as the SM Higgs boson, a CP odd and a pair of charged 
scalars H^. Their masses are given by: 

rnl = —2fil = 2Xiv‘^ 
m]jo = At2 + 

Ttljio — ^2 T XsV 

mjj± = /i 2 + ^Asu^ (3) 

where \l,s are dehned as: 

Al,5 = 2(-^3 + A 4 ± As) (4) 

This model involves 8 independent parameters: hve A, two Hi and v. One parameter is 
eliminated by the minimization condition and the VEV is hxed by the W boson mass. 
Finally, we are left with six independent parameters which we choose as follow : 

{lJ.l,X2,mh,mH±,mHO,mAo} ( 5 ) 
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B. Theoretical and Experimental Constraints 


In order to have vacuum stability, the parameters of the potential need to satisfy the pos¬ 
itivity conditions. Namely, the potential should be bounded from below in all the directions 
of the held space, i.e. should not go to negative inhnity for large held values. We have this 
set of constraints : 

Ai > 0 , A2 > 0 , A3 -|- 2-^/A1A2 > 0 and A3 -|- A4 — [As] > (6) 

We ask that the perturbative unitarity is maintained in variety of scattering processes at 
high energy: scalar-scalar, scalar-vector, vector-vector. Using the equivalence theorem which 
replaces the W and Z bosons by the Goldstone bosons. We hnd a set of four matrices with 
entries are the quartic couplings, the diagonalization of these matrices gives us this set of 
eigenvalues [TB] : 


Cl,2 = A3 ± A4 , 63,4 = A3 ± As 
65,6 = A3 -|- 2A4 ± 3A5 
67,8 = — A2 ± (Ai -|- A2)^ -I- A| 

69,10 = -3Ai - 3 A 2 ± V9(Ai - A 2)2 + ( 2 A 3 + A 4)2 

611,12 = ~Ai — A 2 ± (Ai — A 2 )^ -I- A| (7) 

By requiring that < Stt, we hnd the strongest constraint on Ai ,2 to be: Ai ,2 < Y' 
require that the parameters of the scalar potential remains perturbative, i.e: |Aj| < Stt. 

In order to have an inert vacuum the following constraint should be satished |19j : 

m^±,m^o,m^o > 0 and (8) 

V A1A2 

The extra scalar particles ahects quantum corrections to the W and Z bosons self energies. 
The corrections are parameterized by the oblique parameters, S, T and U [20] which are 
constrained from electroweak precision measurements. Taking the reference Higgs mass as 
nih = 125 GeV and = 173 GeV, the tolerated ranges are found at hxed U=0 m- 

A5 = 0.06 ±0.09, and AT = 0.10 ± 0.07 (9) 
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with correlated factor of +0.91. Where ^nd AT = The 

formulas for AS and AT in the IHDM can be found in Ref (THl 1^ . 

Searches of scalar particles of the IHDM at colliders [23| is not directly performed yet. 
However, several studies [21H26] applied SUSY searches involving two, three or multiple 
leptons with missing transverse energy to the case of IHDM and set some limits on 

the dark Higges. We choose in our analysis : 

> 100 GeV where $ = (10) 

Finally, the magnitude of a possible Higgs boson signal at the LHC is characterized by the 
signal strength modiher, dehned as by : 

^ a{pp —)■ h —)■ Br{h —)■ 77 )™°'^ 

a{pp -)■ hsM -t 77 ) Br{hsM -t 77 ) 

hsM denotes a 125 GeV SM Higgs boson. In our analysis below, while we will show points 
which satisfy theoretical and experimental constraints from our scans, we will highlight the 
points for which Ry^ is consistent with the measured p-yy at the LHC. The latest publicly 
available measurements read [271 EH] 


^ 1.13 + 0.24 (12) 

atlas ^ ('13) 


III. RADIATIVE CORRECTIONS TO TRIPLE HIGGS COUPLINGS hhh AND 

hZZ 

In this section we calculate the one-loop radiative corrections to the trilinear Higgs cou¬ 
pling hhh and hZZ in the SM and IHDM. The correction to those two couplings are part 
of the correction to double Higgs strahlung process e+e“ Zhh which we will discuss in 


^ these dark Higges have the same signature like neutralinos and charginos of the Minimal Supersymetric 
Standard Model (MSSM). 

^ if R or T is DM candidate then its mass could be as low as few GeV. 
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FIG. 1. One loop Feynman diagrams contributing to hhh in the IHDM. S stands for 
and V = W^. 


section V. Those couplings are given at the tree-level by: 


^hhh 


-3ml 


emw 




(14) 

(16) 


As one can see, both couplings hhh and hZZ involve only SM parameters. Those couplings 
receive corrections from one-loop diagrams. The one-loop effects from the SM particles have 
been studied in [301 - 132] for hhh and in [29] for hZZ. These effects are dominated by the top 
quark loops which does not exceed 10% for hhh and 1.5% for hZZ. 

New physics effects to hhh coupling have been analyzed in the context of the Two Higgs 
Doublet Model [30] and the MSSM [32] • It was found that these effects can enhance sig- 
nihcantly this coupling in a wide range of parameter space. Furthermore, these corrections 


depend on the model and hence any deviation from the SM tree level relation (14) by more 
than 10% would be an evidence for the presence of new physics. 

The coupling hZZ have been analyzed in the framework of the two Higgs doublet model 
[3T] and it has been found that the effect is rather small 1% to 2%. 

We have calculated the radiative corrections to the tree level triple Higgs coupling hhh 
and hZZ both in the SM and IHDM in the Feynman gauge including all the particles of the 
model in the loops. The Feynman diagrams from IHDM contributing to \hhh coupling are 
shown in Fig.( [^. 
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The one-loop amplitude are calculated using dimensional regularization. The calculation 
was done with the help of FeynArts and FormCalc [33] packages. Numerical evaluation of 
the one-loop scalar integrals have been done with LoopTools |31|. We have checked both 
numerically and analytically the UV finiteness of the amplitudes. 

In order to do that, we have considered hhh and hZZ at one-loop level: 


i) for hhh, we considered the decay of an off-shell Higgs boson into a Higgs boson pairs 
h*{q) —?• h{ki)h{k 2 ) at the one-loop level. Where q, ki and ^2 cLre the 4-momenta of 
the three particles satisfying on shell conditions ki = k"^ = for hnal state Higgs 
pairs and an off shell condition ^ m\ for the decaying Higgs. 


ii) For hZZ, we follow ref [3T] and write: 


= ml, kl, kl) = (16) 






where ki and k 2 are the momenta of outgoing Z bosons. We assume that the decaying 
Higgs and one of the Z boson are on-shell = ml, kf = m| while the other Z boson is 
off-shell /c 2 = i'^h — "mzY- Using power counting arguments, it is expected that 
receives the highest power contribution of the heavy fermions masses. Therefore, in 
what follow we will take into account only the form-factor to hZZ coupling. 


Since we are dealing with a processes at the one-loop level, a systematic treatment of the 
UV divergences have to be considered. We will use the on-shell renormalization scheme in 
which the input parameters coincide with the physical masses and couplings [35] . In the on 
shell scheme, a redehnition of the helds and parameters is performed. This redehnition cast 
the Lagrangian into a bare Lagrangian and counter-term. The counter-terms are calculated 
by specihc renormalization conditions which allow us to cancel the UV divergences of the 
diagrams with loops. Furthermore, since we have three Higgs as external particles and there 
is no mixing between the SM doublet Hi and the inert doublet H 2 , we do not need to 
renormalize the particle content of the scalar potential of the IHDM. The tree level coupling 


hhh eq. (14) depends only on Higgs mass and the vev as in the SM, then the renormalization 


procedure will be the same as in the SM |35|. We redefine the SM helds and parameters as 





follow: 


my ^ rriy + duly , V = W^,Z 

+ 5ml 

sw —>■ sw + 6 sw 
6 —^ (1 “1“ 
t —)■ t -\- 5t 

^ (l + \5Zzz^ Z^^ + ]^5ZzaA^ 

^ (l + I^Zaa^ + ^SZazZ^ 
h ^ Zl^^h = (^1 + ]^5Zj^ h ( 17 ) 

where sw = sin6*vi/ is the Weinberg angle and t = v{iJi\ — Ain^) is the tadpole which 
is zero at tree level once the minimization condition is used but will receives again hnite 
radiative corrections at the one-loop level. To ensure that the VEV is the same in all orders 
of perturbation theory, it is well known that one need to renormalize the Higgs tadpole: i.e, 
all Higgs tadpole amplitudes T are absorbed into the counter-term 5t. Thus, we put the 
hrst condition: 

f = + T = 0 (18) 

The mass counter-terms are hxed by the on shell conditions [35]: 

ReS^^(m^) = 0 , V = W,Z 

ReZh{m\) = 0 (19) 

The held renormalization constants are hxed by imposing that the residue of the two point 
Green functions to be equal to unity and the mixing y-Z vanish for = m|. While the 
electric charge renormalization constant 5Ze is treated like in quantum electrodynamics and 
is hxed from the e’''e “7 vertex. The renormalized three point function T^+g.^ satisfy at the 
Thomson limit: 

rK+e-.^( A =h = m, = 0) = e 
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Furthermore, the counter-term 6 sw can be obtained from the on-shell dehnition = 1 — 
as a function of Sm^ and 6 mz ■ 

Inserting these redehnitions into the Lagrangian, we hnd the following counter term for 
hhh and hZZ [33] : 


6Chhh — 
6 Chzz = 


—3e^ ml 


2sw m\Y 

SM, = 




6 t 




sw ml 2sw Mwml 

9ci2 _ ^2 ^ 2 




2m^ 


^ + hzn ] 

2 ml, 2 


+ -5Zh + 


hZ^^Z’^ 


( 20 ) 

( 21 ) 


By adding the un-renormalized amplitude for hhh and hZZ to the above corresponding 
counter-terms, one hnd the renormalized amplitudes 


thhh{,q^,rnl,ml) = ,ml,ml) + SChhh 

^hzz{ml, m|, {mh — mzY) = M^h^z + ^ihzz°^ + ^^hzz 

( 22 ) 


which becomes UV hnite. For our numerical illustrations, we dehne the following ratios: 

(23) 

(24) 

Where Thhh is the renormalized vertex. 


AF/i/i/i = 


Re(f^;,;,(g2)) 


^^hZZ = 


A 

-pIDM 
^ hZZ 


hhh 


f SM 
hZZ 


Mi 


IHDM 


Mf 


SM—tree 


Mitzz 


Mf 


SM—tree 


IV. NUMERICAL RESULTS 

A. SM case 

In our numerical analysis, the parameters are chosen as follow : 

mt = 173.5 GeV , mw = 80.3996 GeV , mz = 91.1875 GeV , a = 

127.934 

2 

and the on-shell dehnition of the Weinberg angle: sin^ 6 w = 1 — 

m,z 

In the SM, the dominant contribution to /S.Vhhh{,SM) comes from top quark loops [301132] . 
We have computed the top contribution and shown that it is in perfect agreement with 
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FIG. 2. (Left) ATfihh{SM) as a function of h* momentum q. (Right) ATfizz{SM) as a function 
of A which is the size of deviation from SM triple coupling Xhhh = + A). It is shown: the 

fermionic contribution, the bosonic one as well as the total contribution. 


Ref. [30]. We have also isolated and evaluated the other SM contributions without fermions. 
It turn out that this bosonic contribution is of the order of 5% for large q. 

In Fig. ([^(left), it is illustrated that the top contribution to ATM) is negative before 
the opening of h* —)• ti threshold and also for q > 700 GeV. It is clear that for large q, 
AThhh{SM) is dominated by top-quark contribution. 

In Fig. (|^(right), we show the radiative corrections to hZZ in the SM. We present separately 
the fermionic corrections which are dominated by the top contributions and the bosonic 
contributions. The total corrections to hZZ is of the order of 2%. In this plot, we also shift 
the triple Higgs SM coupling by Af^(l -|- A), where A represent any deviation from 
SM coupling. As one can see from the green line, the sensitivity to A is rather mild. Due 
to custodial symmetry, it is expected that hWW coupling will enjoy similar effect as hZZ 
and that is why we illustrate only the case of hZZ. 


B. IHDM case 

Here, we will show our numerical analysis for the triple coupling of the Higgs in the IHDM 
taking into account: unitarity, perturbativity, false vacuum as well as vacuum stability 
constraints described above. We take the mass of the SM Higgs uih = 125 GeV and the 
masses of the inert particles to be degenerate, i.e : mH± = tuho = = m$. For the other 
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FIG. 3. Left: ATfihh{IHDM) as a function of A 2 . Right: Scatter plot for in the plan 

(m$,/U 2 ) for Q = 300 GeV, Left column show the size of the corrections. and are scanned 
as in eq. (25). 


parameters, we perform the following scan: 

100 GeV < < 500 GeV 

-25 X 10^ (GeV)^ < /i^ < 9 x 10^ (GeV)^ 

4:71 

0 < A2 < y (25) 

We plot in Fig. ([^(left) the relative corrections to the triple coupling as a function of 
A 2 . The theoretical constraints put a limit on A 2 parameter which is A 2 < ^. One can see 
from Fig. [^that the corrections are maximized for A 2 < 2 and decrease for A 2 > 2. In our 
following analysis, We will take A 2 = 2 in order to maximize the effect from A 2 . 

In Fig. (|^(right), we plot the relative corrections to the triple coupling hhh in the plane (m$, 
for a hxed q = 300 GeV and A 2 = 2. One can see that the corrections are very important 
in a large part of the parameter space with an enhancement up to 280% for large values of 
m<i, and negative Furthermore, these corrections are increasing, for a hxed value of m$, 
while fi 2 is decreasing. The maximum of the corrections is reached for ~ —30000 (GeV)^. 
Moderate or very small corrections which can be in the range [—50,50]% are also possible 
for large area of parameter space with low < 300 GeV and any positive It is also 
important to note that LHG constraint from diphoton at the 2cr level exclude light charged 
Higgs 100 < mH± < 175 GeV and negative /i^: left-down corner of the scatter plot. In 
Fig. I^we show the relative corrections AThhhi^HDM) as a function of the momentum of 
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FIG. 4. AThhh{IHDM) as a function of q and where and scanned as in eq. (25). The 

left column shows the values of (left panel) and (right panel) 


the off-shell decaying Higgs q and for and ^2 shown in eq. (25). It is clear that 
for low 100 < < 200 GeV, the corrections are small like in the SM case except at the 

threshold regions where we have ^Thhh ~ —40% dne to the opening of h* —For large 
these corrections could be extremely large exceeding 100 % in large area of parameter 
space. The corrections are amplihed by the opening of the threshold channel h* —)■ <h<I>. 
This is visible on the left panel of Fig. where we can see a kink for q = 1000 GeV which 
correspond to threshold effect with 500 GeV. As it is shown, negative values 

for fil gives large corrections to the triple Higgs coupling. This is because in our assumption 
of taking degenerate Higges 771^0 = m^o = mH± = one can show that A 4 = A 5 = 0, 
-^3 = “ hi) therefore the triple coupling are given by 


hHH = \lv = -(mi — /Xn) 

V 

2 

hAA = \av = -(mi — /in) 

V 

2 


hH+H- = A.n = -{mi 


hi) 


(26) 


It is clear that those couplings gets stronger for negative We now examine the effect 
of the radiative corrections on the triple coupling in the case where the invisible decay 
h —7- HH is open. This is illustrated in Fig. |^left) and right. In Fig. [^left) we impose 
both |Al| < 0.02 required by dark matter constraints as well as best £t limit on the invisible 
branching ratio Br{h —)■ invisible) < 10%(left) and Br{h —>■ invisible) < 20%(right) [36] . 
It is clear from left panel that with dark matter constraint, the size of the corrections and 
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FIG. 5. AF/i/^/i as a function of for q = 300 GeV with /i^ in the range given by the Eq. (25) 
and invisible decay of the Higgs {h HH) is open. Left panel with dark matter constraint which 
restrict \Xl\ to be in the range IAt^I <0.02. Right panel without dark matter constraint. 
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FIG. 6. Left: AThhh{IHDM) as a function of for fixed A 2 = 2 and q = 300 GeV (left), while 
q = 600 GeV in the right. Left column show the range of 


the range of are smaller than in the previous case where the invisible decay was closed. 
The large corrections observed for high mq, are mainly due to the charged Higgs loops. 

In order to exibit the decoupling behavior on the triple Higgs coupling, We increase both 
the range of /i| to be [—10®, 10^] GeV^ as well as the range of G [0.1, 3] TeV. We see that 
the decoupling effect occurs when appropriate combinations of the involved parameters are 
taken large compared to the electroweak scale. We hnd that reaches its maximum 

for mq, ^ 500 GeV and decrease to SM value for large m$. 

In Fig. ([^ we illustrate the IHDM effect on hZZ coupling. Similar to the triple Higgs 


coupling, we £x A 2 = 2 and scan over /i| and mq = mu = = rnH± as in eq. (25). In 

Fig. ([^(left) we show scatter plot for AT^zz in the plan contrarily to the triple 
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FIG. 7. Left: Scatter plot for l^T^zz in the plan with A 2 = 2, left column represent the 

size of the corrections. Right: AT^zz as function of with in the same range as in the left 
panel. Left column show the size of 
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FIG. 8. Feynman diagrams contributing to e^e —)• Zhh at the tree level in SM 


coupling hhh the effect are rather small of the same size as in the SM case. In Fig. ([^ (right) 
we show the corrections to hZZ coupling and also the value of the corresponding within 
2cr range. 


V. RADIATIVE CORRECTIONS TO e+e" ^ hhZ IN THE IHDM 

A. e'’'e“ —7- hhZ in SM 

At LC, the triple Higgs coupling can be probed by double Higgsstrahlung process 
—>■ hhZ depicted in Fig. ([^ and WW fusion process —)• WW*iz9 —)■ hhizV. 

Note that —)■ hhZ process arise in s—channel only and hence its cross section can 

be probed more efficiently at low energies above the threshold (typically between 350 GeV 
and 500 GeV). While, at high energies, for > 700 GeV, the trilinear Higgs-coupling is 
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A 


FIG. 9. Cross section of e"''e —)• Zhh as a function of A where \hhh = ~ 

GeV. 

better probed through the process e'''e“ —?• uPhh assuming the SM. In Fig. (|^ we show the 
Feynman diagrams contributing to e’''e“ —)■ hhZ. When extracting the triple Higgs coupling 
hhh from this process only diagram Fig. ([^(&i) is concerned, the other diagrams 62 , 3,4 are 
considered as a background. 

We illustrate in Fig. ([^ the tree level cross section for e’''e“ —)■ hhZ as a function of A for 
center of mass energy 500 GeV and 800 GeV, where A is the shift of the SM triple Higgs 
coupling Xhhh = + A). It is clear that for A > 0 the cross section is enhanced while 

for A < 0 the cross section is reduced with respect to SM value. 

B. One-loop Corrections to e^e —7- Zhh 

We study in this section the effects of the one-loop radiative corrections to the Higgs tri- 
linear self-coupling calculated in the previous section on the double Higgs-Strahlung process 
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via Z boson exchange. In the context of the SM, the 0{a) electroweak corrections have 
been studied in |37j and it was found that these corrections are of the order 10%. How¬ 
ever, these loop effects can be very large in beyond SM enhancing the total cross section 
by about 2 orders of magnitude in particular in models with extended Higgs sector. As 
outlined above, at the tree level, e+e" —)■ hhZ have four diagrams as depicted in Fig. (|^. 
Only the hrst diagram contribute to the signal while the others are considered as a back¬ 
ground. In our analysis, we include one-loop correction only to the triple Higgs coupling 
hhh which is expected to give sizeable contribution. Only this correction contribute to the 
signal. Therefore, we did not include corrections to the initial state vertex e^e~Z, to the 
self energies of Z-Z and 7 — Z mixing, to the hZZ coupling and also we did not include 
the initial state radiation. In fact these corrections are well known in the SM and are not 
expected to deviate that much in the IHDM as we already show in the previous section for 
hZZ coupling. Moreover, we will not include corrections to e’''e“ —?• hhZ coming from boxes 
and pentagon diagrams. The one-loop amplitude can be written as follow: 


Ad — A^tree T Ad loop 

The squared amplitude at the one-loop level is then : 

\M\^ = + 2Re{M,+ 0(a^) 


(27) 


(28) 


Thus, the cross section is written as: 

+ |Ad|2 (29) 

k=l ^ spin, pola. 

where qi and q2 are the four momentum of the incoming electron and positron, pi,p2 and ps 
are the four momentum of the outgoing particles, and the factor arises from the flux 
of the initial particles. 

For our studies, we dehne the ratio Act by : 


Act 


<^total ^tree 
^tree 


^loop 

^tree 


(30) 


Where cxtotai = o'tree + O'loop- This ratio measure the relative correction of the IHDM to the 
cross section with respect to the tree level result. 
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FIG. 10. (left): Act as a function of for ^/s = 500 GeV values are shown in the right column, 
(right) as a function of A 3 and the right column shows Au in % . 


As stated before, in our analysis we will take into account the theoretical and experimental 
constraints discussed in the second section assuming the parameters to rely in the range given 


by eq. (25). The phase space and evaluation of the one-loop squared amplitude has been 
performed with FormCalc |33] with the help of Loop Tools to evaluate numerically the one- 


loop scalar integrals. We have used the same on-shell renormalization scheme explained in 
the previous section. 


In Fig. 10, we have plotted the ratio Aa versus for center of mass energy y/s = 500 GeV. 
We assume again that: A2 = 2, the dark Higges to be degenerate mriH = mA = 'fnH± = mq, 
and perform a scan over and m$. From this plot one can see that the corrections can reach 
160% for high dark Higgs masses ~ 500 GeV and are negative for low dark Higgs masses 
100 GeV < rriq, < 270 GeV depending on the value of /i^- One can see that the suppression 
of the total cross section can reach —15% for = 140 GeV while the enhancement is 
predominant in most part of the parameter space. 

To understand this, let us remind hrst that in large area of parameter space the correction 
to the triple Higgs coupling /SThhh is positive (see section IV Fig. (|^). Moreover, according 
to the plot Fig. ([^, if this correction is positive this lead to an enhancement of the total 
cross section and vice-versa. 

This explain that in most of the case, the corrections to the cross section are positive and 
conhrm that the behavior of the Aa is consistent with our analysis concerning the trilinear 
Higgs self coupling in the IHDM. We stress that the enhancement of Act is observed in 
a large part of the parameter space and can exceed 100% only in the high mass region 
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> 400 GeV for /i^ < 0. 

We plot in Fig. [Io|( right) the ratio as a function of A3 and showing the relative 
corrections in the left column. For small values of —1 < A3 < 2 (low m$) the corrections 
are quite small. For A3 > 4 (high m$), Act becomes more important and exceed 100%, 
this region corresponds to ~ 0.9 ± 0.02. Note that our results concerning R^^ are in 
agreement with the results of [22] . 


VI. CONCLUSION 

We have computed the radiative corrections to triple Higgs coupling hhh, hZZ coupling 
as well as Zhh in the framework of inert Higgs doublet model taking into account 

theoretical and experimental constraint on the parameter space of the model. The calcula¬ 
tion was done in the Feynman gauge using dimensional regularization in the on-shell scheme. 
In the SM it is known that the top contribution to hhh coupling is of the order 10%, we found 
that the bosonic contribution is somehow signihcant and goes up to 5%. In the IHDM, we 
found that the total radiative corrections to the triple Higgs coupling could be substantial 
exceeding 100% for heavy dark Higgs masses and mH±. We also show that the 

corrections to hhh are decoupling for large m<j, and large In the case of hZZ coupling 
the effect is rather mild and do not exceed 2.5%. We also evaluate radiative corrections to 
the double Higgs strahlung process e+e“ —)■ Zhhhj looking only to the correction to the 
diagram that contribute to the signal i.e the triple coupling hhh. We have shown that the 
correction are also very important. In general, the size of the loop effects, typically large, 
makes their proper inclusion in phenomenological analyses for future e^e“ LC indispensable. 
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APPENDIX 


The Triple Higgs Coupling in the IHDM : Analytical Results 


In this appendix, we present the analytical expression of the Higgs triple coupling at 
the one-loop order with the contribution of the inert scalars only. We use the Feynman 
diagrammatic method. The Feynman diagrams contributing to this process are shown in 
Fig. B By using the dimensional regularization, the amplitude is given by: 


"^l)= ^h±, ^h±) + 


t 2 2 2 2 2 2 \ 

Tia 

(As -|- A4 -|- X5)‘^m\YSw 


X 


Ibevr^ 

2(A3 -|- A4 -|- A5)m^s^ 


Tia 


(As -|- A 4 — X5)‘^rnwSw 


Ibevr^ 

2(As + Xa- A 5 )m^s^ 


X 


ira 


Bo {q^ , , m]jo ) -h 2Ho (m^, m^o, m^o) 

Co(g^ rnl, ml, m|o, m^o, m^o) 
{Bo{q^,m\o,m\a) 2 Ho(m^, m^o, 

^Co{q‘^,ml,ml,m%,m%,m\o) 


( 31 ) 


Where Bq and Co are the Passarino-Veltman functions [38l EH] . 

Following the on-shell renormalization scheme, there are six renormalization constants to 
compute : 6ml,6mly,6m%,6t,6Zf,,6Zh and 6sw SZaa and SZza are the held renormaliza¬ 
tion constant for the photon and Z — 7 mixing respectively and are given by : 




a d 

71 


Hoo(p^m^±,m^±) 


p2=0 


^ZzA 


a{l - 2sly) 
2m|cvi/7r 


(^-Ao{m]j±) + 2Boo{0,mH±,mH 


The electric charge renormalization constant as well as the renormalization constants for the 
W and Z masses are given by: 


= -1(sZaa - —SZza] , ( 33 ) 

/ V Cw / 

Sml^ = (^2Ao (m^±) -F Aq (m^o) + Aq (m^o) 4Hoo (m^, m^±, m^o) ^ (34) 
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5m\ = 


a 


ISvrs^c^ 


-y4o(m^o) - Ao(mjyo) - 2(1 - 2s^) Ao{mH±) + ABoo{mz, m^^o, m^o) 


+4(1 2s^^) Boo^rriz-i Tnjj±, mjj±) 


(35) 


The counter term 5svi/, corresponding to the Weinberg mixing angle, is determined from the 
tree level relation = 1 — Ssw is given by : 


. _ -Cw (Srnj 

- —2 - 2 

2sw V 


(36) 


5ml = 


The counter-terms for the Higgs mass, Higgs held and the tadpole are given by : 

7 ^^ ^2A3Ho(m^±) + (A3 + A4 + A5)Ho(m^o) + (A3 + A4 — A 5 )Ho(m^o)^ 

+ + (A3 + A4 + X 5 fBo{ml,mj^o,mjjo) 

+ (-^3 + A4 — Bo^ml, m\o,m\o)^ 

^ 72 ^? ^^ ^ + A4 + A 5 )^^Ho(p^m^o,m^o) 

9 A 

+ (A3 + A4 — XY^Bq(p‘^, m\o, Tn\o) j 


(37) 


= 


TTLwSw 1 


16e7r^ 327r 


y p2_^2 

^2A3Ho(m^±) + (A3 + A4 + A5)Ho(m^o) + (A3 + A4 — A5)Ho(m^ 


(38) 


(m^orfg) 


The expression for the triple Higgs coupling counter term is given in the third section eq. (20). 
The renormalized triple Higgs coupling is given by : 


-ploop _ -ploop I r-ploop 

^ hhh ~ hhh + hhh 


We have checked that the renormalized amplitude is UV-hnite and furthermore independent 
of the renormalization scale p. We derive an approximate formula for In the limit 
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~ = rnAo = rnH± and find a good agreement with full expression. 




Sanq'^mwsw (2A3 + 2A3A4 + A 4 + A 5 ) ( xi log( + X 2 log(- ^ ^ 2 - 


+ (2A3 + 3A3A4 + 3A3A4 + A4 + 3 (A 3 + A4)A5) rn^s^ log^(-2 ) 


Xi 

2 


X2 


X 1 X 2 


(2A3 + 2A3A4 + A 4 + A 5 ) (ISeTTl^yTTl^ 


s\y{,xi — X 2 ) 

+ dea^m^TT^ - 2m^ + m|)(-l + 2s^) + m5(-2m| + m^(l + (c^ - s^f) 




X ( (-1 + X 2 ) log(—+ (1 - xi) log( ^ ^ 


64m^7rs^(a;i - X 2 ) 


a:ia;2 


X 2 Xi 

aeml{-2m\ + "^^(3 + (cw “ “ 4s^) 


+ 2 + Asi^)) 

X fa;ilog(^^t:^)-a;2log(- 


+ 


X2 

a^m\rn?z'K‘^ 

16em^m|7r^s^ 


Xl 
.2 I ^2 


X ( xf log( -x^ log( 


2mz '^iy(3 + {c\Y '^ly) ^/i( ^ 

I +X2. 


Xl 


X2 


where xi ^2 are given by: 


X 12 = 


1 T ^1 -4m2/g2 


( 40 ) 
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